
global a0
figure(4); % 2D Surface plot

%set(gca,'yLim',[0.1 1]);
%set(gca,'xLim',[0 11]);
set(gca,'XTick', 1:11);
set(gca,'XTickLabel',120:130);
ms=124:0.2:126;
xs=0.5:0.01:0.6;
surf(ms(4:9),xs(1:7),resultsFine(1:7,4:9))
xlabel('Mass point (GeV)');
ylabel('Signal fraction');
zlabel('-ln(L)');

figure(5); % Plot of mass only
hold on;
%set(gca,'yLim',[0.1 1]);
%set(gca,'xLim',[0 11]);
%set(gca,'XTick', 1:11);
%set(gca,'XTickLabel',120:130);
plot(ms,mean(resultsFine(:,:)),'x')
xlabel('Mass point (GeV)');
ylabel('-ln(L)');
ms=124:0.2:126;
a0=polyfit(ms(4:9),mean(resultsFine(:,4:9)),2);
xpoints=124:0.01:126;
plot(xpoints,quadraticFunc(xpoints));

figure(6); % Plot of ratio only
hold on;
%set(gca,'yLim',[0.1 1]);
%set(gca,'xLim',[0 11]);
%set(gca,'XTick', 1:11);
%set(gca,'XTickLabel',120:130);
xs=0.5:0.01:0.6;
plot(xs(1:6),resultsFine(1:6,5),'x')
xlabel('Signal ratio');
ylabel('-ln(L)');

a0=polyfit(xs(3:8),resultsFine(3:8,5)',2);
xpoints=0.5:0.001:0.55;
plot(xpoints,quadraticFunc(xpoints));